The c2d Spitzer spectroscopic survey of ices around low-mass 
young stellar objects. IV. NH3 and CH3OH 
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ABSTRACT 

NH3 and CH3OH are key molecules in astrochemical networks leading to the formation of more 
complex N- and 0-bearing molecules, such as CH3CN and CH3OCH3. Despite a number of recent 
studies, little is known about their abundances in the solid state. This is particularly the case 
for low-mass protostars, for which only the launch of the Spitzer Space Telescope has permitted 
high sensitivity observations of the ices around these objects. In this work, we investigate the 
~ 8 — 10 fj,ia region in the Spitzer IRS (InfraRed Spectrograph) spectra of 41 low-mass young 
stellar objects (YSOs). These data are part of a survey of interstellar ices in a sample of low-mass 
YSOs studied in earlier papers in this series. We used both an empirical and a local continuum 
method to correct for the contribution from the 10 /im silicate absorption in the recorded spectra. 
In addition, we conducted a systematic laboratory study of NH3- and CHsOH-containing ices to 
help interpret the astronomical spectra. We clearly detect a feature at ~9 /im in 24 low-mass 
YSOs. Within the uncertainty in continuum determination, we identify this feature with the 
NII3 1^2 umbrella mode, and derive abundances with respect to water between ~2 and 15%. 
Simultaneously, we also revisited the case of CH3OH ice by studying the 1/4 C-0 stretch mode 
of this molecule at ^9.7 /Ltm in 16 objects, yielding abundances consistent with those derived 
bv iBoogert et al.l ( 2008[) (hereafter paper I) based on a simultaneous 9.75 and 3.53 /im data 
analysis. Our study indicates that NH3 is present primarily in Il20-rich ices, but that in some 
cases, such ices are insufficient to explain the observed narrow FWHM. The laboratory data point 
to CII3OII being in an almost pure methanol ice, or mixed mainly with CO or CO2, consistent 
with its formation through hydrogenation on grains. Finally, we use our derived NH3 abundances 
in combination with previously published abundances of other solid N-bearing species to find that 
up to 10-20 % of nitrogen is locked up in known ices. 

Subject headings: infrared: ISM — ISM: molecules — ISM: abundances — stars: formation — astro- 
chemistry 
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1. Introduction 

Ammonia and methanol are among the most 
ubiquitous and abundant (after H2 and CO) 
molecules in space. Gaseous NH3 and CH3OH 
are found in a variety of environments such 
as infrared dark clouds, dense gas surrounding 
ultra-compact H II regions, massive hot cores, 
hot corinos, and comets. Solid CH3OH has 
been observe d in the ices surrounding massive 



YSOs 



[e^ 



Schutte et al.l Il991t iDartois et al 



19991: iGibb et al.ll2004l) and more recently towar d 



low-mass protostars (Pontop pidan et al.l 120031) 
The presence of solid NH3 has been c l aimed 
toward massive YSOs onlv (ILacy et al. 1998t 



Dartois et al.ll2002l:lGibb et al.ll2004HGiirtler et ~ 



20021 ) ■ with the exception of a possible detec- 
tion in the low-ma ss object IRAS 03445-1-3242 
(JGiirtler et al.l 120021 ). However, the se detections 



are st ill controversial and ambiguous (jTaban et al 
20031) . 



Both molecules are key participants in gas-grain 
chemical networks resulting in the formation of 
more complex N- and 0-bear ing molecules, such as 



CH3CN and CH3OCH3 (e.g. iRodeers fc CharnlevI 
20011 ). Moreover, UV processing of NH3- and 
CH30H-containing ices has been proposed as a 
way to produce amino- a cids and other complex 
organic molecules (e .g., iMunoz Caro fc Schutte 



2003tlBernstein et al.ll2002l: lOberg et al.ll2009[) . In 

addition, the amount of NH3 in ices has a direct 
impact on the content of ions such as NHj and 
OCN~ , which form reactive intermediates in solid- 
state chemical networks. A better knowledge of 
the NH3 and CH3OH content in interstellar ices 
will thus help to constrain chemical models and 
to gain a better understanding of the formation of 
more complex, prebiotic, molecules. 

During the pre-stellar phase, NH3 is known to 
freeze out on grains (if the c ore remains starless 
long enough - iLee et al.ll2004 ). Moreover, CH3OH 
is known to have gas-phase abundances with re- 
spect to H2 in hot cores/corinos that are much 
larger than in cold dense clouds: ■^ (1 — 10) x 
10~^ vs < 10^^, with the former values most 
likely repr esenting evaporated ices in warm re- 
gions (e.g. iGenzel et al.l Il982l: iBlake et aLlll987 : 



Federman et al.lll99d) . Together, these findings 
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suggest that ices are an important reservoir of 
NH3 and CH3OH and that prominent features 
should be seen in the absorption spectra toward 
high- and low-mass protostars. Unfortunately, as 
summarized in Table [TJ NH3 and CH3OH bands, 
with the exception of the 3.53 /im CH3OH feature, 
are often blended with deep water and/or silicate 
absorptions, complicating unambiguous identifica- 
tions and column density measurements. This is 
particularly true for NH3 whose abundance deter- 
mination based on the presence of an ammonium 
hydr ate feature at 3.47 /xm remain s controversial 
(e.g. iDartois fc d'Hendecourtll200l[ ). Nonetheless, 
it is important to use all available constraints to 
accurately determine the abundances of these two 
molecules. Despite the overlap with the 10 /im 
silicate (Si-O stretch) feature, the NH3 V2 um- 
brella mode at ~9 /xm (^^1110 cm~^) offers a 
strong intrinsic absorption cross section and ap- 
pears as the most promising feature to determine 
the abundance of this species in the solid phase. 
Moreover, the CH3OH v^ C-O stretch at ~9.7 /xm 
(~1030 cm~^) provides a good check on the valid- 
ity of the different methods we will use to subtract 
the 10-/im silicate absorption, since the abundance 
of this molecule has been accurately determined 
previously from both the 3.53 and 9.75 /um fea- 
tures (see Paper I). 

More detailed spectroscopic information is par- 
ticularly interesting for low-mass protostars as 
the ice composition refiects the conditions dur- 
ing the formation of Sun-like stars. Such de- 
tections have only become possible with Spitzer, 
whose sensitivity is necessary to observe low lu- 
minosity objects even in the nearest star-forming 
clouds. The spectral resolu tion of the Spitzer In- 
frared Spectrograph (IRS, iHouck et al] |2004| ) of 
AA/A ~ 100 in this wavelength range is compa- 
rable to that of the Infrared Space Observatory 
(ISO) PHOT-S instrument but lower than that of 
the ISO-SWS and other instruments used to iden- 
tify solid NH3 toward high luminosity sources. 
The spectral appearance of ice absorption fea- 
tures, such as band shape, band position and in- 
tegrated band strength, is rather sensitive to the 
molecular environment. Thus, the interpretation 
of the astronomical spectra should be supported 
by a systematic laboratory study of interstellar ice 
analogues containing NH3 and CH3OH. Changes 
in the lattice geometry and physical conditions 



Table 1 
Selected near- and mid-infrared features of NH3 and CH3OH. 



Mode 



A (/im) V (cm ^) Problem 



NH3 features: 






z/3 N-H stretch 


2.96 


3375 


1/4 H-N-H bend 


6.16 


1624 


IV2 umbrella 


9.00 


1110 


CH3OH features: 






V2 C-H stretch 


3.53 


2827 


vq & f3 -CH3 deformation 


6.85 


1460 


V7 -CH3 rock 


8.87 


1128 


1/4 C— stretch 


9.75 


1026 


Torsion 


14.39 


695 



Blended with H2O (0-H stretch, 
3.05 Aim/3275 cm^i) 
Blended with H2O (H-O-H bend, 
5.99 Atm/1670 cm"!), HCOOH 
Blended with silicate 



Blended (e.g. with NH;^) 
Weak; blended with silicate 
Blended with silicate 
Blended with H2O libration mode 



Note. — The bold-faced lines indicate the features studied here. 



Note.— The nomenclature for the NH3 and CH3OH vibrational modes are adopted from 
Herzberej (|l945l ). 



of an ice are directly reflected by variations in 
these spectral properties. In the laboratory, it is 
possible to record dependencies over a wide range 
of astrophysically relevant parameters, most ob- 
viously ice composition, mixing ratios, and tem- 
perature. Such laboratory data exist for pure 
and some H oO-rich NH3- and CH.-^QH-con t aining 



ices (e.g.. Id'Hen decour t fc Allamandqla 



Hudgins et al.lll99 3: Kcrk hof et al.lll999l : lTaban et all 
20031 ). but a systematic study and comparison 
with observational spectra is lacking. 

In principle, the molecular environment also 
provides information on the formation pathway of 
the molecule. For example, NH3 ice is expected to 
form simultaneously with II2O and CH4 ice in the 
early, low-density molecular clo ud phase from hy- 



droge nation of N atoms (e.g., iTielens fc Hagen 
19821 ). In contrast, solid CII3OII is thought 
to result primarily from hydrogenation of solid 
CO, a process which has been confirmed to be 
rapid at low temper atures in several labora- 
tory experinients (e.g. IWatanabe fc Kouchil 12002 ; 



Hidaka et ai1l2004t iFuchs et al.ll2009l ). A separate. 
water-poor layer of CO ice is often found on top of 
the water-rich ice layer in low-mass star-forming 
regions due to the 'catastrophi c' freeze-out of gas- 



phase CO at high dens ities (jPontoppidan et al 



2OO3I : IPontoppidanI I2OO6I ) . Hydrogenation of this 
CO layer s hould lead to a nearly pure CII3OH ice 
layer (e.g.. ICuppen et al.l 12009). which will have 
a different spectroscopic signature from that of 
CH3OH embedded in a water-rich matrix. The 
latter signature would be expected if CH3OH 
ice were formed by hydrogenation of CO in a 
water-rich environment or by photoprocessing of 
H2 0: C0 ice mixtures, ano ther proposed route 
(e.g., iMoore fc Hudsonlll998l) . 

Here, we present Spitzer spectra between 5 and 
35 /zm of ices surrounding 41 low-mass protostars, 
focusing on the ~ 8 — 10 /im region that contains 
the 1/2 umbrella and 1^4 C-0 stretch modes of NH3 
and CH3OH, respectively. Th is work is the fourth 
paper in a series of ice studies (IBoogert et al. l2008t 
Pontoppidan et al.l l2008t lOberg et al.l l2008j car- 
ried out in the context of the Spitzer Legacy Pro- 
gram "From M olecular Cores to Planet- Forming 
Disks" ("c2d": IEws et al.|[2003l) . In Section [H 
we carry out the analysis of the Spitzer data in 
8 — 10 /im range. In Section[3l we present the lab- 
oratory data specifically obtained to help interpret 



the data that are discussed in Section ID Finally, 
we conclude in Section [5] with a short discussion 
of the joint astronomy-laboratory work (including 
the overall continuum determination). 

2. Astronomical observations and analysis 

The source sample consists of 41 low-mass 
YSOs that were selected based on the presence 
of ice absorption features. The entire sample 
spans a wide range of spectral indices a = —0.25 
to +2.70, with a defined as dlog(AFx)/dlog(A), 
where d indicates the derivative, and F\ repre- 
sents all the photometric fluxes available between 
A = 2.17 ^m (2MASS Kg-band) and A = 24 ^m 
{Spitzer/MIPS band). In the infrared broad-band 
classification scheme, 35 out of 41 objects fall in 
the embedded Class O/I category {a > 0.3). The 
remaining 6 obje cts are flat-spectrur a type objects 
(-0.3 <a< 0.3: lGreeneet al.lll994l ). 
Spitzer/IRS spectra (5-35 /xm) were obtained as 
part of the c2d Legacy program (PIDs 172 and 
179), as well as a dedicated open time program 
(PID 20604), and several previo usly published 
GTO spectra (jWatson et a l.||2004|). We refer the 
reader to Table 1 and Section 3 of iBoogert et al 



(|20Q8I ) for the source coordinates and a description 
of the data reduction process (including overall 
continuum determination) . 

As mentioned previously, spectral signatures in 
the ~ 8 — 10 /im region are dominated by the Si- 
O stretching mode of silicates. The overall shape 
as well as the sub-structure of the silicate feature 
depend on grain size, mineralogy, level of crys- 
tallinity. These effects are degenerate and so these 
different factors cannot be easily separated. For 
example, large grains and the presence of SiC both 
produ ce a shoulder at 11.2 /im (e.g 



Min et al 



20071 ). Therefore, trying to fit the 10 fim silicate 
feature by determining the composition and size 
of the grains is a complex process. For this rea- 
son, we use two alternative methods to model the 
silicate profile and extract the NH3 (and CH3OH) 
feature(s) from the underlying silicate absorption. 

2.1. Local continuum 

The first method uses a local continuum to fit 
the shape of the silicate absorption. For this, we 
flt a fourth order polynomial over the wavelength 



ranges 8.25-8.75, 9.23-9.37, and 9.98-10.4 /zm, 
avoiding the positions where NH3 and CH3OH ab- 
sorb around 9 and 9.7 fim. These fits are shown 
with thick black hnes in Fig. [T] After subtrac- 
tion of the local continuum from the observations, 
we fit a Gaussian to the remaining NH3 and/or 
CH3OH feature, when present, as shown in Fig. [2] 
The results of the Gaussian fits are listed in Ta- 
bic lA.ll of Appendix El 

2.2. Template 

The second method assumes that the 8-10 /nm 
continuum can be represented by a template sil- 
icate absorption feature, selected among the ob- 
served sources. A comparison of the results ob- 
tained using a template to those obtained using 
a simple local continuum provides an estimate of 
the influence of the continuum choice on the shape 
and depth of the NH3 and CH3OH features. The 
templates were chosen using an empirical method. 
Upon examination of the 10 /xm feature of the en- 
tire sample, the sources could be separated into 
three general categories, depending on the shape of 
the wing of the silicate absorption between '^8 and 
8.7 /xm (which we will refer to as the 8 /xm wing): 
(i) sources with a straight 8 ^m wing (Fig. [3]-a) , 
(ii) sources with a curved 8 //m wing (Fig. [3]-b) , 
and (iii) sources with a rising 8 /xm wing ("emis- 
sion" sources, Fig. 131-c) . 

Note that, since radiative transfer in the 8- 
10 /im region can be complicated by the presence 
of silicate emission, we only consider sources that 
are the least affected by emission, that is those 
falling in one of the first two categories. Never- 
theless, non-rising silicate profiles might still suffer 
from the presence of emission. To try and estimate 
the impact of this potential effect, we used two sil- 
icate emission sources from iKessler-Silacci et al.l 
(|2006h . and subtracted these emission profiles 
from our absorption profiles, assuming that the 
emission represented 10 to 50% of the observed 
absorption. After removal of a local continuum, 
we determined the integrated optical depths of the 
NH3 and GH3OH features in the spectra corrected 
for emission, and compared these to the integrated 
optical depths of the uncorrected spectra. We find 
that the difference can be up to a factor of two and 
therefore identify this possible presence of under- 
lying emission as the largest source of uncertainty 
in our abundance determinations. 



For each of the straight and curved 8 /im 
wings, two sources (in order to test for template- 
dependent effects) were selected as possible tem- 
plates for the silicate feature. The selection crite- 
ria were: (i) a silicate feature as deep as possible 
to minimize the effects of silicate emission and (ii) 
little NH3 and CH3OH signal, as estimated after 
subtraction of a local continuum. Additionally, we 
added to this l i st the GCS3 spectrum observed by 
Kemper et al.l (|2004l ) toward the Galactic Center. 
The spectra of these templates in the 8-10 ^m 
region are displayed in Fig. |4l 

For all the other sources in our sample, the best 
template was determined by scaling the possible 
templates to the observed optical depth at differ- 
ent wavelengths (8.75, 9.30, 9.37, 9.70, 9.98 ^J,m) 
and finding the combination (template-fscaling 
point) that gave the least residuals over the same 
wavelength ranges used to estimate the local con- 
tinuum (8.25-8.75, 9.23-9.37, 9.98-10.4 ^m). The 
result of this process is displayed for each source 
in the top part of Fig. [1] where the best template 
is shown by a grey line. The bottom panels of 
Fig. [I] show sources for which no reasonable tem- 
plate could be found, as well as emission sources, 
in which case only the local continuum is overlaid. 
As in the case of the local continuum method, the 
spectra obtained after subtraction of the templates 
are shown in Fig.[2j Taken together, NH3 features 
are detected in 24 out of 41 sources. 

The top panel of Figure [5] shows that the 
CH3OH feature is not affected by the continuum 
choice, whereas the width of the NH3 band is 
somewhat sensitive to this choice, especially if 
there is no CH3OH absorption, in which case the 
local continuum yields a wider NH3 profile. For 
both continua, there is a clear feature around 9 
fim, which we attribute to NH3, with the char- 
acteristics and limitations given and discussed in 
the following sections. 

2.3. NH3 ice column densities and abun- 
dances 

Gaussian fits were performed to the NH3 
and/or CH3OH features when present, and de- 
rived parameters for NH3 are listed in Table lA.ll 
(Appendix |X| . Table [2] gives the column densi- 
ties derived for NH3 for each of the two meth- 
ods employed to determine the continuum, us- 
ing a band strength of 1.3x10^^^ cm molecule^ ^ 



for the NH3 i/2 u mbrella mode appropriate for 
a wate r-rich ice (Id'Hendecourt fc Allamandolal 
[l986 : .Kerkhof et al.lll999l ). The two methods gen- 
erally agree to within a factor of 2 or better. A 
similar factor of <2 overall micertainty is esti- 
mated for those sources for which only the local 
continuum has been used. 

The position of the NH3 1^2 umbrella mode is 
very close to that of the v^ CHa-rock mode of 
CH3OH. As illustrated by our laboratory data (see 
Section [3]) , sources with an absorption depth at 
~9.7 /im (CO-stretch mode of CH3OH) at least 
twice as large as the absorption depth at '-~^9 /im 
(blend of CH3-rock mode of CH3OH and NH3 
umbrella mode) have a significant contribution to 
the 9-/im integrated optical depth from the CH3- 
rock mode of CH3OH. In these cases (sources fol- 
lowed by an asterisk in Table [2] and in Table lA.ll 
of Appendix |Al , we performed the following cor- 
rection: we scaled a H20:CH30H=9:1 labora- 
tory spectrum to the observed optical depth of 
the CO-strctch mode of CH3OH, determined the 
integrated optical depth of the CH3-rock mode 
of CH3OH in that scaled spectrum, and sub- 
tracted it from the total observed optical depth 
at 9 /xm. This correction is justified by the fact 
that the H2O:CH3OH:NH3=10:4:l spectrum, a 
typical interstellar abundance mixture, is well re- 
produced around 8-10 /im by a combination of 
H20:CH30H=9:1 and H20:NH3=9:1 (see SecEJ. 

The inferred NH3 ice abundances range from 
< 1% to 15% with respect to H2O ice, excluding 
the abnormally high value of EC 82. When consid- 
ering all values (except that of EC 82) determined 
with the local continuum method, this relative 
abundance is centered on 5.3% with a standard 
deviation of 2.0%. If we use values determined 
with the template method whenever available, we 
find a mean of 7.0±3.2%. Either way, within the 
errors, this i s sini ilar to what was obtained by 
Qberg et al.l (12008") for CH4 (4.7±1.6%), another 



ice component that should form via hydrogena- 
tion. For 6 out of the 8 sources where both NH3 
and CH4 are detected, the NH3-to-CH4 abun- 
dance ratio is slightly larger than 1 ('-^1.2). Based 
on elemental abundance ratios, one would expect 
NH3/CII4 smaller than 1, but since two thirds of 
the carbon is in refractory grains and some frac- 
tion of the gaseous CO locked up in CO at the 
ice formation threshold, NH3-to-CIl4 ratios larger 



than 1 are consistent with both NII3 and CH4 
being formed by hydrogenation of N and C, re- 
spectively. Further comparison between these and 
other ice species will be addressed in an upcoming 
paper (Oberg et al. in preparation). 

Regarding CII3OII, we only report values for 
the Gaussian parameters and derived column den- 
sities in the appendix (see Table [5l2|) . to show that 
the numbers we obtain in this independent study 
are consistent with those reported in Paper I. Our 
recommended abundances are those from paper I, 
based on the combined 9.75 and 3.53 /im analy- 
sis. The inferred CH3OH abundances range from 
< 1% to > 25% with respect to H2O ice, indi- 
cating significant CH3OH/NH3 abundance varia- 
tions from source-to-source. Such relative abun- 
dance variations can already be clearly seen from 
the changing relative depths of the 9.0 and 9.7 /im 
features (see also Paper I). Thus, NII3 and CH3OH 
ice are likely formed through different formation 
pathways and/or in different ice environments. 
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Fig. 1. — (Top) Local continuum (thick blue/black lines) and template (red/grey lines) fits to all sources for which a template 
could be found. (See Section 12.21 for details) — [Bottom) Local continuum fits to emission sources or sources for which no 
reasonable template could be found. 
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Fig. 2. — (Top) Residual after removal of local continuum and template fits for all sources for which a template could be 
found. (See Section 12.21 for details) — (Bottom) Residual aftetiremoval of local continuum fits for emission sources or sources 
for which no reasonable template could be found. 
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Fig. 3. — Examples illustrating the three shapes of the 8 jxth wing shown by the thick grey line: (a) straight, 
(b) curved, and (c) rising. 
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Fig. 4. — Silicate features of the sources used as templates for a straight 8 /im wing (left), curved 8-/im wing 
(middle), and GCS3 (right). The bottom panels of each plot are the residuals after removal of the local 
continuum shown in grey in the top panels. The optical depth scale is kept fixed for comparison. These 
sources are selected to have no or at most weak NH3 and CH3OH absorptions. 



Table 2 
nh3 column densities" and abundances with respect to h2o ice^ 



Source 



NH3, local 



% HaO" 



NH3, template 



xlO'^ 



% HsG" 



Template 



Scaling point 
/iin 



IRAS 03235+3004 


6.83 ( 0.98) 


4.71 ( 1.00) 


8.94 ( 1.03) 


6.17 ( 1.20) 


IRAS 12553 


9.30 


L1455 IRS3 


0.57 ( 0.23) 


6.21 ( 3.51) 


1.41 ( 0.27) 


15.37 ( 6.86) 


GCS3 


9.37 


IRAS 03254+3050 


2.44 ( 0.39) 


6.66 ( 1.37) 


4.58 ( 0.49) 


12.52 ( 2.10) 


IRAS 12553 


10.40 


Bl-b* 


~7.3 


~4.2 


~9.8 


~5.6 


IRAS 12553 


9.70 


IRAS 04108+2803 


1.23 ( 0.24) 


4.29 ( 1.03) 


2.07 ( 0.39) 


7.21 ( 1.69) 


IRAS 23238 


9.70 


HH 300 


0.90 ( 0.22) 


3.46 ( 0.90) 


2.23 ( 0.37) 


8.60 ( 1.65) 


DC Tau B 


9.70 


IRAS 08242-5050 


4.77 ( 0.46) 


6.13 ( 0.85) 


4.41 ( 0.54) 


5.66 ( 0.89) 


IRAS 12553 


9.70 


IRAS 15398-3359 


8.73 ( 1.18) 


5.90 ( 1.77) 


13.80 ( 1.35) 


9.33 ( 2.65) 


IRAS 12553 


9.70 


B59 YS05 


4.92 ( 0.72) 


3.53 ( 0.88) 


6.37 ( 0.99) 


4.57 ( 1.17) 


CrA IRS7 A 


9.70 


2MASSJ17112317-272431 


13.10 ( 1.06) 


6.70 ( 0.54) 


20.60 ( 2.76) 


10.58 ( 1.42) 


IRAS 23238 


9.70 


SVS 4-5* 


~2.4 


~4.3 


~5.8 


~10.3 


GCS3 


8.75 


R CrA IRS 5 


0.91 ( 0.23) 


2.54 ( 0.67) 


1.49 ( 0.31) 


4.15 ( 0.92) 


IRAS 12553 


9.70 


RNO 15"= 


0.80 ( 0.21) 


11.58 ( 3.18) 


- 


- 


- 


- 


IRAS 03271+3013 


4.90 ( 0.88) 


6.37 ( 1.86) 


- 


- 


- 


- 


Bl-a 


3.46 ( 0.69) 


3.33 ( 0.98) 


- 


- 


- 


- 


L1489 IRS 


2.31 ( 0.30) 


5.42 ( 0.96) 


- 


- 


- 


- 


IRAS 13546-3941 


0.94 ( 0.16) 


4.56 ( 0.87) 


- 


- 


- 


- 


RNO 91 


2.03 ( 0.30) 


4.78 ( 0.81) 


- 


- 


- 


- 


IRAS 17081-2721 


0.86 ( 0.16) 


6.54 ( 1.39) 


- 


- 


- 


- 


EC 74"= 


1.00 ( 0.29) 


9.35 ( 3.13) 


- 


- 


- 


- 


EC 82 


1.22 ( 0.14) 


31.31 ( 6.65) 


- 


- 


- 


- 


EC 90 


0.67 ( 0.20) 


3.94 ( 1.24) 


- 


- 


- 


- 


EC 92* 


~0.5 


~3.0 


- 


- 


- 


- 


CrA IRS 7 B* 


~3.0 


~2.8 


- 


- 


- 


- 


L1014 IRS 


3.72 ( 0.91) 


5.20 ( 1.43) 


- 


- 


- 


- 


CK4 


0.84 ( 0.13) 


5.37 ( 0.86) 


- 


- 


- 


- 



3-(7 upper limits 



LDN 1448 IRSl 


0.20 


4.15 


IRAS 03245+3002 


17.28 


4.40 


L1455 SMMl 


15.10 


8.29 


IRAS 03301+3111 


0.24 


5.93 


Bl-c 


11.93 


4.04 


IRAS 03439+3233 


0.31 


3.10 


IRAS 03445+3242 


0.47 


2.09 


DG Tau B 


0.47 


2.05 


IRAS 12553-7651 


0.61 


2.04 


Elias 29 


0.28 


0.93 


CRBR 2422.8-342 


0.52 


1.23 


HH 100 IRS 


0.46 


1.89 


CrA IRS7 A 


0.97 


0.89 


CrA IRAS32 


5.44 


10.35 


IRAS 23238+7401 


1.60 


1.24 



Note. — Sources in bold were used as templates. Uncertainties quoted in parenthesis arc statistical errors from the Gaussian fits 
while absolute errors are up to a factor of 2. 

^Derived using a bandstrength of 1.3xl0~ cm molecule" . 

Using the H2O ice column densities listed in Paper I. 

'^Values are likely upper limits (see Scction |4.2| for details). 

*Sources with T9 7.^111 > 2 X T9 ounii for which an estimated contribution from the CHs-rock mode of CH3OH was subtracted (see 
text for details). 
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3. Laboratory work and analysis 

The band profiles presented in Fig. [5] contain 
information on the ice environment in which NH3 
and CH3OH are located, and thus their formation 
and processing history. To extract this informa- 
tion, a systematic laboratory study of the NH3 
and CH3OH features in a variety of ices has been 
carried out. Specifically, three features between 8 
and 10 /im have been analyzed: 

1. the NH3 V2 umbrella mode, at '^9.35 /im or 
1070 cm^^ in pure NH3 ice, and with band 
=1.7xl0~^^ cm molecule^"'^ 



Wavelength (/^m) 
9.0 9.5 



strength Annrp 

(|d'Hendecourt fc Allamandolalll986l) . 



the CH3OH Vii CO- stretching mode, at 
--9.74 /im or 1027 cm^^ in pure CH3OH ice. 



and with Ap„rR — 1.8x10' 



^^ cm molecule 



(|d'Hendecourt fc Allamandolal[l986l) 



the CH3OH vi CH3 rocking mode, at 
-8.87 /xm or 1128 cm^^ in pure CH3OH ice. 



and with A. 



"'"•R , 



(|Hudgins et al.l ll993) 



1.8x10 ° cm molecule 



It should be noted that, as mentioned in the 
above list, the quoted positions are for pure ices 
only and therefore slightly deviate from the astro- 
nomical values given in Table [TJ 

This laboratory study targeted pure, binary 
and tertiary interstellar ice analogs consisting of 
different mixtures of H2O, NH3, CH3OH, CO 
and CO2, the major ice components. All mea- 
surements were performed under high vacuum 
conditions (~ 10~^ mba r) using an ex perimen- 
tal approach des c ribed in lGerakines et al. (1995) . 



Bouwman et al.l (|2007n and lOberg et al.l (|2007n 



The ice spectra were recorded in transmission 
using a Fourier transform infrared spectrome- 
ter covering 25-2.5 /im (400-4000 cm^^) with 
1 cm~^ resolution and by sampling relatively 
thick ices, typically several thousands Langmuiiu 
(L) thick. These ices were grown at a speed of 
~10^® molecules cm~^ 
controlled Csl window. 



on a temperature- 



^One Langmuir corresponds to a pressure of 10~^ torr for 
1 second and measures the exposure of a surface to ad- 
sorption of gases. One Langmuir is equivalent to about 
10^^ molecules cm~'^. 



o o ; 

S " 0.020 r 



H20;CH30H:NH3=10:4:1 

H20:CH30H-9:1 

H20:NH3=9:1 

H20:CH30H-9: 1 +H20;NH3-9: 1 




1150 1100 1050 

Wavenumber (cm"') 



Fig. 5. — Example of a reduced laboratory spec- 
trum (solid black line) for a H20:CH30H:NH3= 
10:4:1 ice mixture at 15 K, in the 8-10 /xm / 960- 
1220 cm~^ range. This spectrum can be approx- 
imated as the sum (solid green/dark grey line) of 
H2 0:CH30H=9:1 (solid red/light grey line) and 
H20:NH3=9:1 (dash-dot blue/grey line) . The bot- 
tom plot is the difference between the two, show- 
ing that the feature at 9 /im (blend of NH3 and 
CH3OH CHs-rock modes) is well reproduced by 
the sum of the two individual signatures. This fig- 
ure also illustrates the fact that the positions of 
the features in mixed ices differ from that in pure 
ices (see list at the beginning of this section). 



A typical reduced spectrum for an ice mixture 
containing H20:CH30H:NH3= 10:4:1 at 15 K is 
shown in Fig. [S] Since band profiles and strengths 
change with ice composition and also with tem- 
perature, the three fundamentals mentioned above 
were investigated as a function of temperature 
ranging from 15 to 140 K with regular tempera- 
ture steps for a number of binary and tertiary mix- 
tures (listed in Appendix |B| . An IDL routine was 
used to determine the location of the band max- 
imum, FWHM and integrated absorbance of the 
individual absorption bands. For the asymmet- 
ric NH3 vi umbrella mode the band position has 
been determined by the maximum absorbance and 
for the symmetric profiles the spectral parameters 
have been determined from Gaussian fits of base- 
line subtracted spectra. The resulting absolute 
frequency uncertainty is of the order of 1 cm~^. 
The measurements are presented in Table IB. 31 of 
Appendix B, and are included in the Leiden labo- 



11 



Wavelength (jum) 
B.5 9.0 i 



10.0 10.5 
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-0.5 
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Fig. 6.— {Left) FTIR ice spectra of the ^nhs mode for pure NH3, a H20:NH3=1:1 and a HaOrNHg^g:! 
mixture at a temperature of 15 K. At the low frequency side of the spectrum the H2O libration mode 
(centered around 770 cm^^, or 13 /im) starts showing up for the H20-containing mixtures. — (Right) 
Temperature effect on a H20:NH3=9:1 mixture: decreasing FWHM with increasing temperature. 



ratory databasqj. 

NH3 and CH3OH both have the ability to form 
hydrogen bonds in water-rich matrices, so it is not 
surprising that the band profile changes compared 
with pure ices bec ause of the various molecular 
interactions (e.g., Id'Hendecourt fc Allamandola 



19861) 

In addition to profiles, band strengths can 
change with environment and with temperature, 
as discussed f or the cases of CO and CO2 in water- 

(199i); 



rich ices in 



(|20071) : iBouwman et al.l 



Kerkhof et al.| . 
(I2OO7I) 



Oberg et al 



Figure |6] shows 
how the NH3 V2 umbrella mode absorption max- 
imum shifts from 1070 cm~^ (9.35 /im) for pure 
NH3 ice to 1118 cm^^ (8.94 ^m) for an astro- 
nomically more realistic H20:NH3=9:1 (hereafter 
9:1) mixture, for which the FWHM and integrated 
band strength also change significantly. For exam- 
ple, the band strength is lowered in the 9:1 mixture 
to 70% of its initial value in pure NH3 ice. This 
is in good ag reement with previou s experiments 
performed bv iKerkhof et al.l ( 1999I ). The spectral 
appearance also depends on temperature; for the 
9:1 mixture a temperature increase from 15 to 
120 K results in a redshift of the peak position 
from 1118 to 1112 cm'^ (8.94 to 8.99 ^m) and 
the FWHM decreases from 62 to 52 cm^^ (0.50 to 
0.42 /im) (see Fig.[7l). The NH3 bandstrength, on 
the other hand, does not show any temperature 



'^http : //www ■ strw . leldenunlv . nl/~lab/ 1 



dependence. 

If NH3 is in a water-poor environment with CO 
and/or CO2, the z^2 peak position shifts to the 
red compared with pure NH3, to as much as 1062 
cm-i (9.41 ^m). The FWHM is not much affected 
whereas the band strength is lowered by 20%. Be- 
cause of the intrinsically large difference in band 
maximum position between NH3 in a water-poor 
and water-rich environment, the astronomical ob- 
servations can distinguish between these two sce- 
narios. 

Methanol-containing ices have been studied in 
a similar way (see Fig. |8]). The weakly absorbing 
vr CH3 rocking mode at ~1125 cm~^ (8.89 ^m) 
is rather insensitive to H2O mixing, but the v^ 
CO stretch vibration shifts to the red from 1028 
to 1020 cm~^ (9.73 to 9.80 /im) when changing 
from a pure CH3OH ice to a H20:CH30H=9:1 
mixture. In the latter spectrum the CH3OH 1/4 
CO stretch mode needs to be fitted with a dou- 
ble Gaussians. A substructure appears for a tem- 
perature of 80 K (right panel of Fig. |S]) while for 
even higher temperatures, a clearly double peaked 
structure bec omes visible (as pr eviously seen in 
e.g. Fig. 2 of lSchutte eraLlll99l[ ). This splitting 
hints at different physical sites and has been pre- 
viously ascribed to type I I clathrate formation in 
the ice (JBlake et al.lll99ll) . 

When CH3OH is mixed with CO, the band 
maximum shifts from 1028 to 1034 cm"! (9.73 
to 9.67 /im) when going from a 9:1 to a 1:9 
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Fig. 7. — A plot indicating the changes in peak position (left) and FWHM (right) of the NH3 V2 umbrella 
mode as a function of temperature in a 9:1 H20:NH3 ice. 
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Fig. 8. — [Left) Spectra of the CH3OH i^co modes and t^cHs modes for pure CH3OH, a H20:CH30H=1:1, 
a H20:CH30H=9:1 and a C0:CH30H=1:1 ice mixture at a temperature of 15 K. — {Right) Temperature 
effect on the CO-stretch mode of a H20:CH30H=9:1 mixture. 



CH30H:C0 mixture. When 50% or more CO 
is mixed in, the CH3OH v^ CO stretch mode 
starts to show a shoulder and cannot be fitted cor- 
rectly by a single Gaussian component (see Fig. [9]). 
Such a two-component profile would not be recog- 
nized, however, at the spectral resolution and sig- 
nal/noise of our Spitzer data, so for the compari- 



son between laboratory and observational data a 
single Gaussian is used. Overall, the shifts of the 
CH3OH ly'i mode between water-rich and CO-rich 
mixtures are much smaller than in the case of the 
NH3 1/2 mode. 

The effect of CH3OH on the 4.7^m Vi stretch 
mode of CO has also been investigated. The band 
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Fig. 9. — Spectra of CHsOHiCO mixtures in the 
range of the methanol CO stretch mode and the 
methanol CH3 rock mode. A small blue shift to- 
gether with a clear substructure are seen upon 
mixing in more CO. 

maximum shifts from 2139 cm~^ (4.68 /im) for 
the nearly pure 9:1 C0:CH30H mixture to 2136 
and 2135 cm~^ for the 1:1 and 1:9 mixtures, re- 
spectively. The CO band located at 2136 cm~^ is 
often referred to as CO residing in a polar, mainly 
H2O ice, environment. Clearly, the polar CH3OH 
molecules can also contribute to CO absorption 
at 2136 cm~^ when intimately mixed in an astro- 
nomical ice. 

Binary mixtures of NH3 and CH3OH have been 
studied as well. The CH3OH modes behave very 
much as they do in a pure methanol ice, but 
the NH3 1^2 umbrella mode is clearly suppressed. 
Its integrated absorbance is readily reduced to 
70% of the integrated absorbance of pure NH3 
in a CH30H:NH3 = 1:1 mixture and becomes even 
lower for a 4:1 binary composition. The NH3 
band also broadens compared to pure NH3 or 
H20:NH3 mixtures and strongly overlaps with the 
CO stretching mode of CH3OH, to the level that 
it becomes difficult to measure. 

A qualitative comparison with the astronom- 
ical data (see Section |4]) indicates that neither 
pure NH3, CH3OH, nor mixed CH30H:NH3 or 
H2 0-diluted binary ices can simultaneously ex- 
plain the different NH3 profiles in the recorded 
Spitzer spectra. Thus, a series of tertiary mix- 
tures with H20:CH30H:NH3 in ratios 10:4:1, 
10:1:1 and 10:0.25:1 have been measured, be- 
cause CH3OH is the next major ice compo- 



Fig. 10. — Normalized spectra of the 

CH3OH t/4 C-0 mode (right panel), 

and NH3 z^2 umbrella mode (left panel) 
for a H2O:CH3OH:NH3=10:0.25:l, 

a H2O:CH3OH:NH3=10:l:l and a 

H2O:CH3OH:NH3=10:4:l mixture at a tem- 
perature of 15 K. These mixture ratios span the 
range of observed interstellar column density ra- 
tios. Spectra were normalized to better show the 
changes in band maximum position and FWHM of 
each feature. Spectra of a H20:CH30H==9:1 and a 
H20:NH3=1:1 mixture were offset and overlaid in 
light grey in the right and left panel, respectively. 
In the case of H2O:CH3OH:NH3=10:4:l, the NH3 
1/2 umbrella mode is heavily blended with the 
CH3OH 1/7 CH3 rocking mode, so that the dark 
grey line actually shows the Gaussian fit to the 
underlying NH3 feature, whereas the full 9-/im 
feature is shown in black. 



nent. These ratios roughly span the range of 
observed interstellar column density ratios. In 
Fig. Uni the spectra of H20:CH30H:NH3 ter- 
tiary mixtures are plotted and compared to bi- 
nary H20:CH30H and H20:NH3 data. The NH3 
1/2 umbrella mode shifts slightly to the blue in 
the presence of both H2O and CH3OH, with an 
absorption maximum at 1125 cm~^ (8.90 /im) for 
the 10:4:1 H20:CH30H:NH3 mixture (compared 
to 1118 cm-i (8.94 /im) in the H20:NH3=9:1 mix- 
ture). The peak intensity of the NH3 1/2 umbrella 
mode band in this tertiary mixture is small com- 
pared with that of the CII3OII CH3 rock mode, 
but its integrated intensity is a factor of two larger 
because of the larger NH3 width. 

The 1/4 C-0 stretching vibration profile of 
CH3OH in the tertiary mixture does not differ 
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Fig. 11. — FWHM and band maximum position of the 
NH3 feature measured in the laboratory mixtures at 15 K 
("Lab.", top panel) and in the 5piteer spectra ("Astro.", bot- 
tom panel). In the top panel, the symbols are labelled with 
their corresponding mixtures ; for the NH3:CH30H mixtures 
(orange triangles), an increasing symbol size is indicative of 
increasing CH3OH content. In the bottom panel, open and 
filled squares indicate values obtained with the template and 
local continuum method, respectively. In both panels, the 
dash-dot polygons delimitate the parameter space of FWHM 
and positions corresponding to H20-rich mixtures. 
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Fig. 12. — Same as Fig.[TT]but for CH3OH. 



much from the binary values for the highest water 
content. The position of the absorption maxi- 
mum is also only marginally affected by the tem- 
perature. The FWHM decreases from 30 cm^^ 
(0.29 iim) for the 10:4:1 mixture to 22 cm^^ 
(0.21 ^m) for the 10:0.25:1 mixture. 

Besides H2O, other species may also be re- 
garded as potential candidates for changing the 
spectral appearance of the NH3 and/ or CH3OH 



featu res. Chemically linked is HCOOH (Bisschop et al. 



20071 ) which unfortunately cannot be deposited in 
the present setup because of its reactive behavior 
when mixed with NH3. Tertiary mixtures with 
CO and CO2, two other important constituents 
in interstellar ices, have been measured (see Ap- 
pendix |B]) but here the differences are small com- 
pared with the observed binary water-rich or CO- 



rich mixtures, and do not offer an alternative ex- 
planation. 

4. Comparison between astronomical and 
laboratory data 

4.1. 8-10 fim range 

The FWHM and band positions of the NH3 and 
CH3OH features measured in the laboratory and 
astronomical spectra are shown in Figs. [11] (for 
NH3) andim (for CH3OH). For the YSOs, the 
values obtained after removal of the silicate ab- 
sorption (see Section [2|) using the local continuum 
method are indicated by filled squares, whereas 
those obtained from the template method are plot- 
ted with open squares. Note that the presence 
of significant amounts of CH3OH may artificially 
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lower the inferred NH3 V2 width in CH3OH rich 
sources (indicated with * in Table[2|) because of the 
contribution of the narrower vj CHs-rock mode. 

Regardless of the method used to subtract 
the continuum, or the type of source (CH3OH- 
rich/poor), we find that the observational band 
positions and FWHM of the V2 NH3 umbrella 
mode absorptions vary, within the errors, between 
8.9 and 9.1 /im and between ~0.2 and 0.5 /im, re- 
spectively. These position and width are not well 
simultaneously reproduced by any of the investi- 
gated mixtures. Regarding the positions, those 
measured in water-rich ice mixtures are the clos- 
est, whereas the positions in pure NH3 or CO/CO2 
rich ices are too far away to be representative of 
the astronomical positions. The derived Spitzer 
FWHM values range between 0.23 and 0.32 ^m 
(except for Bl-b : 0.39 /im), when using the lo- 
cal continuum method, not depending on whether 
the target is CH30H-rich or -poor. For the tem- 
plate method, CHsOH-rich sources generally tend 
to have a narrower inferred FWHM, 0.3-0.5 /J,m, 
contrary to what would be expected if the NH3 
mode is contaminated by the CHs-rock feature. In 
any case, most of these widths are still narrower 
than the laboratory FWHM values. To investigate 
further the effect of the continuum on the posi- 
tions and widths of the bands, we performed the 
following alternative analysis to check whether a 
continuum could be found that would yield NH3 
and CH3OH features with parameters within the 
laboratory measurements. To do that, we fitted 
the data between 8.25 and 10.4 fim with a function 
that is the sum of a 4'^ order polynomial and two 
gaussians ; positions and widths of the gaussians 
were constrained with limits taken from the labo- 
ratory data of binary water mixtures (8.9-8.95 ^m 
for the NH3 position, 0.42-0.52 /im for its width; 
9.67-9.77 /im for the CH3OH position, 0.2-0.3 /im 
for its width). As illustrated in Fig. [TH we found 
that the continuum derived in this way is different 
from those determined via the other two methods. 
This result supports the fact that the difference 
between astronomical and laboratory data could 
be attributed to the uncertainty in the continuum 
determination. 

Taking the above considerations into account. 
Figs . [m and [T2l suggest that the template method 
for subtraction of the 10 /im silicate absorption 
is more consistent with the laboratory measure- 



ments, but both methods probably miss some 
weak NH3 absorption features in the broad line 
wings where they blend with the continuum at 
the S/N of the data. If so, the too small line 
widths inferred from the data (most probably due 
to the uncertainty in the continuum determina- 
tion) would mean that we have underestimated 
NH3 abundances by a up to a factor of 2. 

The observational band position and FWHM 
of the CH3OH features derived with either the lo- 
cal continuum or template method are clustered 
around 9.7-9.75 /im, with the exception of R 
CrA IRS 5 at 9.66 /im. Similarly the FWHM 
of the CH3OH features are all very similar be- 
tween -0.22 and 0.32 /^m, except for R CrA IRS 
5 with 0.39 /im. These values agree (with a few 
exceptions) with the values obtained from the lab- 
oratory spectra. Note that the observed posi- 
tions of the CH3OH feature are all on the low 
side of the laboratory range. Since the position 
of this feature shifts to higher wavelengths with 
increasing water content, the observed low values 
could therefore indicate that CH3OH and H2O are 
not well mixed and that there exists a separate 
CH3 OH-ri ch component, as suggested in previous 



work (e.g. Ponto ppidan et al.ll2003t ISkinner et al 
19921 ). Alternatively, the low values could be due 



to the presence of CO as indicated by the CH3OH 
feature shift to 9.70 /im in CH30H:C0=1:1. Both 
interpretations would be consistent with the bulk 
of the CH3OH formation coming from hydrogena- 
tion of a CO-rich layer, rather than photochem- 
istry in a water-rich matrix. However, the shift 
from the water-rich mixtures is small, and some 
water-rich fraction cannot be excluded with the 
current spectral resolution. 

4.2. The 3 and 6 /im ranges 

iDartois fc d'HendecourtI (|200l[ ) discussed the 



possibility of a 3.47 /im absorption band which 
could be related to the formation of an am- 
monia hydrate in the ice mantles: they found 
that if this band were mostly due to this hy- 
drate, then ammonia abundances would be at 
most 5% with respect to water ice. Considering 
the fact that our derived abundances are larger 
than 10% in some sources, it is necessary to in- 
vestigate the effect of such a high abundance on 
the ammonia features in other spectral ranges. 
For this, depending on the NHs-to-CHsOH abun- 
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Fig. 13. — Example of the effect of contin- 
uum determination. Upper panel: Spectrum of 
IRAS 03235+3004 in the 8.25-10.4 /im region over- 
laid with the overall fit (smooth purple line - sum of 
a 4**^ order polynomial and two gaussians with posi- 
tions and widths constrained by the values measured 
in the laboratory spectra of H20:NH3 and CH3OH 
mixtures), the 4*'^ order polynomial obtained by the 
overall fit (red) and the local continuum previously 
determined (green). The purple histogram is the dif- 
ference between the data and the overall fit. — Bot- 
tom panel: Residuals after subtraction of the two 
continua, in respective colors. The yellow and blue 
lines are H20:NH3=9:1 laboratory spectra scaled to 
the red and green residuals, respectively. Note the 
good agreement between the feature extracted with 
the red continuum and the lab data (orange) show- 
ing that astronomical and laboratory data of NH3 
ice mixtures can be consistent if a slightly different 
continuum determination is adopted. 



dance ratio observed in the Spitzer spectra, we 
scale one of the following laboratory spectra to the 
9 ^ira NH3 feature: H20:NH3=9:1, H20:NH3=4:1, 
H2O:CH3OH:NH3=10:l:l,H2O:CH3OH:NH3=10:4:l. 
Figure [14] illustrates the comparison between the 
Spitzer and scaled laboratory spectra for the rel- 
evant wavelength ranges for a couple of sources, 
while Fig. lC.lf a and -b (see Appendix IBJ) show the 
comparison for all sources where NH3 was tenta- 
tively detected. 

For further comparison, we also overplotted in 
Fig. [TJ] and IC.ll the following spectra: (i) the 
pure H2O ice spectrum d erived from th e H2O 
column density quoted in iBoogert et al.l (J2008I ) 
(deep blue); and (ii) for sources with 3 /im data, 
the pure H2O spectrum scaled to the optical depth 
of the 3-/im feature of the mixed ice laboratory 



spectrum (purple-dotted). The difference between 
this scaled pure water spectrum and the mixed ice 
spectrum gives an indication of the contribution 
of ammonia features around 3.47 and 6.1 /im. 

We then determined the contributions from the 
NH3 features to the integrated optical depths of 
the 3 and 6 /im bands and to the optical depth of 
component C2, a feature at 6.0 — 6.4 fiiri arising 
from a blend of several species, including NH3, 
H2O, CO2, HCOO~ (see Paper I for more details). 
These contributions are reported in Table [3l 

Figure [H ED and Table [3] show that (i) the 
scaled laboratory spectra generally do not overes- 
timate the observed absorption features, and (ii) 
for most sources, the presence of NH3 at the level 
we determine from the 9 /im feature does not ex- 
plain by itself the depth of the C2 component and 
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Fig. 14. — Comparison of astronomical data (VLT or Keck measurements at short wavelengths, IRS Spitzer 
observations elsewhere) and laboratory spectra in selected wavelength ranges: 2.0 — 4.5 ^m (left panels), 
5.2 — 7.5 /im (middle panels) and 8.2 — 10.2 /im (right panels, silicate absorption subtracted via the template 
method). Error bars are indicated in the bottom-right corner. Overlaid in red and green are laboratory 
spectra corresponding to H2O:CH3OH:NH3=10:4:l and H20:NH3=9;1, respectively, scaled to the 9-/im 
NH3 umbrella mode, and smoothed to the Spitzer resolution in the ~ 5 — 10 ^m range. The dark blue line 
represents the pure water laboratory spectrum scaled to the water column density taken in paper I. The 
dotted purple line corresponds to a pure water spectrum scaled to the 3-/im water feature of the mixed ice 
spectrum, showing the contribution of NH3 features around 3.47 and 6.1 /im. Finally, the red dashed line in 
the right panel of SVS 4-5 represents a H20:CH30H=9:1 laboratory spectrum scaled to the 9.7-/im CH3OH 
CO-stretch mode: this gives an indication of the contribution of the 9-/im CH3OH CH3-rock mode to the 
total 9-/im feature. The laboratory spectra are recorded at 15 K unless indicated differently. 
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Table 3 
nh3 contribution to the 3 and 6 ^m bands 



Source 


/THaO.a.o 

f-^3.0 


f 


iix,3.0 
T3.0 


rl785 

Jl562 ^H20 
.1785 
J1562 


rl785 . 
J 1562 '^^^ 
rl785 ^ 
Jl562 


/tnh, 

(•1785 
Jl562 


3,6.16 
THsO 


''"NH3,6.16 


TC2 


Sources with template 


IRAS 03235+3004 


- 




- 




0.50 




0.24 




0.02 


0.61 


IRAS 03254+3050 


0.73 




1.30 




0.56 




0.92 




0.12 


1.72 


IRAS 04108+2803 


0.70 




0.67 




0.58 




0.53 




0.06 


0.49 


HH 300 


0.70 




0.57 




0.50 




0.39 




0.05 


0.45 


IRAS 08242-5050 


0.76 




0.72 




0.50 




0.45 




0.06 


0.46 


IRAS 08242-5050 


0.76 




0.56 




0.50 




0.35 




0.05 


0.36 


2MASSJ17112317-272431 


- 




- 




0.69 




0.53 




0.05 


4.23 


SVS 4-5 


0.91 




0.94 




0.42 




0.29 




0.00 


0.08 


R CrA IRS 5 


0.85 




0.42 




0.63 




0.29 




0.03 


0.21 




Sources with no 


associated template 










RNO 15 


0.80 




1.97 




0.53 




1.23 




0.16 


0.45 


IRAS 03271+3013 


- 




- 




0.36 




0.44 




0.05 


0.60 


Bl-a 


- 




- 




0.67 




0.43 




0.03 


0.57 


L1489 IRS 


0.78 




0.88 




0.60 




0.56 




0.04 


0.83 


RNO 91 


0.94 




0.94 




0.53 




0.45 




0.04 


0.53 


IRAS 17081-2721 


0.65 




0.95 




0.62 




0.75 




0.05 


1.64 


EC 74 


0.95 




2.34 




0.57 




1.18 




0.09 


0.76 


EC 92 


0.90 




0.35 




0.38 




0.10 




0.00 


0.01 


CrA IRS 7 B 


- 




- 




0.81 




0.19 




0.00 


0.08 


L1014 IRS 


- 




- 




0.62 




0.55 




0.06 


0.34 



Note. — A dash indicates that the ratio was not calculated due to the high noise in the 3-/j.m spectrum. 
Parameters are: 

/rHgOjS.O — integrated optical depth of pure water at 3 /im, determined from the column density of paper I and a 
band strength of 2.0x10"^^ cm^^. 

/''"3.O1 / '''mix, 3.0 — integrated optical depth over the entire 3 /xm region for, respectively, the considered source and 
the corresponding laboratory mixture (selected from the NH3 feature at 9 /im). 



// 



1785 , 
562 



THoO 



n 



1785 
562 



II 



1785 
562 



integrated optical depth of, respectively, pure water, source spectrum, and labo 



ratory mixture, between 1562 and 1785 cm~ (5.6 to 6.4 /im). 

/ '''NH3,6.i6, ''"NH3,6.i6 — integrated and peak optical depth of the 6.16 /im feature of ammonia obtained afte 

subtraction of a pure water spectrum scaled to the optical depth at 3 /im of the laboratory mixture. 

tq2 — peak optical depth of the C2 component from paper I. 
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of the red wing of the 3 fim band. Hence, our 
inferred NH3 abundances up to 15% from the 9.7 
/Ltm data are not in conflict with the lack of other 
NH3 features. The only exceptions are two sources 
(RNO 15 and EC 74), for which the scaled mixed 
ice spectrum exceeds the data in the 3-/j,m range. 
In the case of RNO 15, the NH3 abundance could 
have been overestimated due to the contribution of 
the CH3OH CH3-rock feature at -9 ^im. For EC 
74, this overestimate and the presence of emission 
weakens the identification of the NH3 signature. 
In both cases, the quoted NH3 abundances should 
be considered as upper limits. 

Overall, our reported NH3 abundances are up 
to a fact or of three larger than the upper limits de- 
rived by JDartois^&^^d^Hendecouo (|200^ Firstly, 
let's recall that the conclusions in their study and 
in ours are drawn from the analysis of different 
samples. Secondly, Dartois & d'Hendecourt made 
an assumption that does not apply to our sample: 
indeed, they considered a grain size distribution 
including also scattering from larger grains, pro- 
ducing an enhanced 3-/im wing, whereas the re- 
sults presented here can be taken as representative 
of NH3 absorption from small grains. It is beyond 
the scope of this paper to investigate the effects of 
grain size di stribution and scattering in as much 
detail as did iDartois &: d'HendecourtI (|200l[ ) . 



• Identification 
OCN- (e.g. 



of substantial amounts of 



van Broekhuizen et al. 200 



4.3. Nitrogen ice inventory 

The confirmation of the presence of relatively 
large amounts of solid NH3, up to 15%, in inter- 
stellar ices solves a long-standing problem. In- 
deed, the detection of solid NII3 has remained 
elusive and/or controversial, despite a number of 
clues suggesting its presence: 

• High cosmic abundance of atomic nitrogen : 
iWiVn = 7.76 X IQ-^ (Savage & SembacJ 



19961), only a factor of a few below those 
of oxygen and carbon. Here Nh indi- 
cates the total number of hydrogen nuclei, 

NH^N{R)+2N(ti2)- 

• High abundances of gaseous NH3 of ^nh^ /./V"h-> ^ 
IQ-^-lO"'^ in the O rion-KL nebula (JBarrett et al 
I977I: iGenzel et all [l982i) and in other hot 
cores such as G9.62-F0.19, G29.96 -0.02, 
G31.41-H0.31 ("Cesaroni et al.' '1994*), and 



2005 ) and NHJ" in ices (e.g . Schuttc & Khanni 



2003t iBoogert et aLllJOoi) : considering that 



3 

umal 



these ions form via reactions involving NH3 , 
the non-detection of solid NH3 would be 
puzzling. 

Our results can be used to draw up a possible 
nitrogen budget. Assumin g A^Hoo/A^h ^ 5 x 10~^ 



(jPontoppidan et al.l 120041 : iBoogert et aD l2004l ). 



and average abundances w.r.t. H2O of 5.5% for 
NH3 (see glSl), 7% for NH + (from Table 3 of Paper 



I), a nd 0.6% for OCN (jvan Broekhuizen et al 



G10. 47-H0.03 ( Cesaroni et al.iil994iOsorio et al. 
20091 ). 



20051 ). then the NH3, NH| and OCN" abundances 
with respect to total H are 2.8, 3.5, and 0.3 xlO~^ 
respectively. This corresponds to, respectively, 
3.4, 4.4 and 0.4% of the atomic nitrogen cosmic 
abundance so that, in total, about 10% of the cos- 
mically available nitrogen would be locked up in 
ices, leaving solid and gaseous N2, N and HCN 
as other substantial nitrogen carriers. The main 
uncertainty in this determination is the adopted 
H2O ice abundance with respect to total H; in sev- 
eral sources this may well be a factor of 2 larger, 
leading to about 20% of the nitrogen accounted 
for in ices. 

5. Conclusion 

We have analyzed in detail the 8-10 fim range 
of the spectra of 41 low-ma ss YSOs obtain e d wit h 
Spitzer and presented in iBoogert et al.l (|2008l ). 
The sources are categorized into three types: 
straight, curved and rising 8 /Ltm silicate wings, 
and for each category template sources with lit- 
tle or no absorption from ices around 9-10 /xm 
have been determined. This has led to two ways 
of subtracting the contribution from the 10 /im 
silicate absorption: first, by determining a local 
continuum, and second, by scaling the templates 
to the optical depth at 9.7 /xm. The two methods 
give consistent band positions of the NH3 features, 
but the resulting widths can be up to a factor of 
two larger using the template continuum method. 
Taking into account the uncertainty in continuum 
determination, NH3 ice is most likely detected 
in 24 of the 41 sources with abundances of ^2 
to 15 % w.r.t. H2O, with an average abundance 
of 5.5±2.0 %. These abundances have estimated 
uncertainties up to a factor of two and are not 
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inconsistent with other features in the 3 and 6 /itm 
ranges. CH3OH is often detected as well, but the 
NH3/CH3OH abundance ratio changes strongly 
from source to source. Our inferred CH3OH col- 
umn densities are consistent with the values de- 
rived in paper I. 

Targeted laboratory experiments have been car- 
ried out to characterize the NH3 and CH3OH pro- 
files (position, FWHM, integrated absorbance). 
Comparison with the observational data shows 
reasonable agreement (within ~1%) for the posi- 
tion of the NH3 feature in H2 0-rich ices, but the 
observed widths arc systematically smaller than 
the laboratory ones for nearly all sources. The 
silicate template continuum method gives widths 
that come closest to the laboratory values. This 
difference in width (i.e. widths derived from as- 
tronomical spectra smaller than those in the labo- 
ratory spectra) suggests that the NH3 abundances 
determined here may be on the low side. 

The CH3OH profile is most consistent with a 
significant fraction of the CH3OH in a relatively 
pure or CO-rich phase, consistent with its forma- 
tion by the hydrogenation of CO ice. In contrast, 
the most likely formation route of NH3 ice remains 
hydrogenation of atomic N together with water 
ice formation in a relatively low density molecular 
phase. Finally, the nitrogen budget indicates that 
up to 10 to 20 % of nitrogen is locked up in known 
ices. 
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A. Parameters of Gaussian fits 
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Table A.l 
Parameters of Gaussian fits to the NH3 feature. 



Source 



NH3, local 



A (pm) FWHM (^m) 



NH3, template 



A (pm) FWHM (/^m) 



RNO 15 

IRAS 03271+3013 

Bl-a 

L1489 IRS 

IRAS 13546-3941 

RNO 91 

IRAS 17081-2721 

EC 74 

EC 82 

EC 90 

EC 92* 

CK4 

CiA IRS7 B* 

L1014 IRS 



9.05±0.02 
8.96±0.02 
8.98±0.02 
9.02±0.01 
8.99±0.02 
8.98±0.01 
8.97±0.02 
9. Olio. 02 
8.94±0.01 
8.95±0.02 
8.99±0.02 
8.99±0.02 
9.04±0.01 
9.03±0.02 



0.20±0.04 
0.25±0.04 
0.25±0.04 
0.27±0.03 
0.27±0.03 
0.26±0.03 
0.24±0.04 
0.22±0.05 
0.29±0.03 
0.18±0.05 
0.25±0.05 
0.32±0.04 
0.26±0.03 
0.25±0.05 



0.04±0.01 
0.20±0.02 
0.14±0.02 
0.09±0.01 
0.03±0.00 
0.08±0.01 
0.04±0.00 
0.05±0.01 
0.04±0.00 
0.04±0.01 
0.03±0.00 
0.03±0.00 
0.15±0.01 
0.15±0.02 



Tpcak 



IRAS 03235-1-3004 8.93±0.02 0.28±0.03 0.23±0.02 8.93±0.01 0.29±0.03 0.30±0.02 

L1455 IRS3 8.99±0.03 0.24±0.07 0.02±0.01 9.02±0.02 0.38±0.05 0.04±0.01 

IRAS 032544-3050 9.04±0.01 0.25±0.03 0.10±0.01 8.99±0.01 0.38±0.03 0.12±0.01 

Bl-b* 9.05±0.03 0.39±0.06 0.25±0.02 9.07±0.03 0.40±0.06 0.31±0.03 

IRAS 04108+2803 8.99±0.02 0.25+0.04 0.05+0.01 9.05+0.03 0.47+0.06 0.04+0.01 

HH 300 9.01+0.02 0.23+0.05 0.04+0.01 9.06+0.02 0.45+0.06 0.05+0.01 

IRAS 08242-5050 9.02+0.01 0.31+0.03 0.15+0.01 9.05+0.01 0.30+0.03 0.15+0.01 

IRAS 15398-3359 8.96+0.01 0.29+0.03 0.30+0.02 8.98+0.01 0.33+0.03 0.41+0.02 

B59 YS05 8.95+0.01 0.27+0.03 0.18+0.02 8.89+0.02 0.34+0.04 0.18+0.02 

2MASSJ17112317-272431 8.99+0.01 0.30+0.02 0.43+0.02 9.02+0.02 0.50+0.05 0.41 + 0.04 

SVS 4-5* 9.00+0.01 0.26+0.03 0.16+0.02 9.01+0.01 0.30+0.03 0.26 + 0.02 

R CrA IRS 5 9.05+0.02 0.21+0.04 0.04+0.01 9.00+0.03 0.36+0.06 0.04+0.01 



Note. — Uncertainties arc statistical errors from the Gaussian fits. 

*Sources with T9.7^rn > 2 X rg.o^im, for which the contribution from the CH3OH CHs-rock mode is significant. Since 
the latter and the NH3 umbrella mode were difficult to disentangle, a single fit was performed (the reported parameters) 
and the integrated optical depth of the ammonia feature was then obtained from the total integrated optical depth at 
9 fimhy subtracting the estimated contribution of the CH3OH CH3-rock mode (sec i|2.2t . 
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Table A. 2 

Parameters of Gaussian fits to the CH3OH C-O stretch mode (after subtraction of the continuum with the local 

and/or template method), and CH3OH column densities (or 3-0- upper limits). 



to 



Source 




Local continuum 






Template 


continuum 




Paper I 




A 


FWHM 


'''peak 


X 


A 


FWHM 


"^pcak 


X 


X 




(^m) 


(/jin) 




(% H2O) 


(/^m) 


(/^m) 




(% H2O) 


(% H2O) 


IRAS 03235+3004 


9.74+0.02 


0.26+0.03 


0.35+0.04 


4.40+1.04 


9.74+0.02 


0.25+0.04 


0.31 + 0.04 


3.84+0.99 


4.20+1.20 


L1455 IRS3 


9.78+0.01 


0.14+0.03 


0.03+0.01 


3.67+1.80 


9.78+0.02 


0.26+0.04 


0.04+0.01 


7.71+3.46 


<12.5 


IRAS 03254+3050 








< 5.4 








< 5.4 


< 4.6 


Bl-b 


9.71+0.01 


0.30+0.03 


1.19+0.11 


14.15+3.16 


9.71+0.01 


0.28+0.03 


1.21+0.11 


13.75+3.12 


11.20+0.70 


IRAS 04108+2803 








< 2.7 


9.75+0.00 


0.06+0.04 


0.04+0.03 


0.58+0.62 


< 3.5 


HH 300 








< 4.7 


9.74+0.00 


0.19+0.12 


0.01+0.01 


0.78+0.52 


< 6.7 


IRAS 08242-5050 


9.70+0.01 


0.27+0.03 


0.25+0.02 


6.12+1.01 


9.70+0.01 


0.29+0.03 


0.24+0.02 


6.39+1.09 


5.50+0.30 


IRAS 15398-3359 


9.73+0.01 


0.28+0.03 


0.77+0.06 


10.26+3.02 


9.73+0.01 


0.30+0.03 


0.75 + 0.06 


10.69+3.14 


10.30+0.80 


B59 YS05 








< 1.2 








< 1.2 


< 1.3 


2MASSJ17112317-272431 


9.75+0.02 


0.23+0.04 


0.13+0.02 


1.03+0.22 








< 2.0 


< 3.2 


SVS 4-5 


9.74+0.01 


0.28+0.03 


0.77+0.06 


26.38+6.17 


9.74+0.01 


0.31 + 0.02 


0.83 + 0.06 


31.50+7.12 


25.20+3.50 


R CrA IRS 5 


9.66+0.01 


0.39+0.03 


0.07+0.00 


5.68+0.60 


9.66+0.02 


0.39+0.04 


0.07+0.00 


5.51+0.72 


6.60+1.60 


RNO 15 


9.65+0.03 


0.44+0.07 


0.02+0.00 


11.13+2.16 


















< 5.0 


IRAS 03271+3013 








< 4.3 


















< 5.6 


Bl-a 








< 2.4 


















< 1.9 


L1489 IRS 


9.78+0.02 


0.10+0.03 


0.03+0.01 


0.44+0.22 


















4.90+1.50 


IRAS 13546-3941 








< 2.0 


















< 3.9 


RNO 91 


9.77+0.01 


0.11+0.03 


0.05+0.01 


0.87+0.32 


















< 5.6 


IRAS 17081-2721 








< 6.6 


















3.30+0.80 


EC 74 








<13.5 


















< 9.3 


EC 82 








<24.6 


















<14.2 


EC 90 


9.70+0.01 


0.32+0.03 


0.05+0.00 


6.91+0.99 


















6.80+1.60 


EC 92 

CK4 

CrA IRS7 B 


9.73+0.01 


0.30+0.02 


0.09+0.01 


11.16+1.46 


















11.70+3.50 


9.70+0.01 


0.33+0.02 


0.36+0.02 


7.74+1.56 


















6.80+0.30 


L1014 IRS 


9.69+0.03 


0.38+0.08 


0.10+0.01 


3.61+0.99 


















3.10+0.80 



Note. — This tabic shows that CH3OH column densities obtained in this paper arc consistent with those in Paper I, which arc our recommended values. 
Note. — Uncertainties arc statistical errors from the Gaussian fits. 



B. Additional information on laboratory data 
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Table B.3 
Ice composition, band maximum position ("peak position"), FWHM and band strength 
relative to the pure ice (^ircl.), listed for a set of ice mixtures under investigation. 





Ice 


mixture 






Molecule 


Peak 


position 


FWHM 


Arol. 




Mode 


NH3 


CH30H 


H2O 


CO 


C02 


cin^^ 


^m 


cm^"^ 


/^m 




















NH3 


1070 


9.341 


66 


0.577 


1 


1^2 


umbrella 







0.11 








NH3 


1076 


9.291 


70 


0.605 


1 


V-2 


umbrella 







1 








NH3 


1100 


9.091 


77 


0.637 


1 


VI 


umbrella 







9 








NH3 


1118 


8.947 


62 


0.496 


0.7 


VI 


umbrella 







10 


1 





NH3 


1124 


8.897 


53 


0.420 


0.7 


vi 


umbrella 







1 


1 





NH3 


1094 


9.144 


75 


0.627 


1 


V2 


umbrella 







10 





2 


NH3 


1122 


8.916 


57 


0.453 


0.8 


Vl 


umbrella 







1 





1 


NH3 


1098 


9.108 


82 


0.681 


0.9 


V2 


umbrella 










1 


1 


NH3 


1062 


9.414 


66 


0.586 


0.8 


V2 


umbrella 




4 











NH3 


1129" 


8.856° 


108° 


0.849 


0.4° 


V2 


umbrella 




4 











CH3OH 


1029 


9.722 


30 


0.283 


- 


Ui 


C-O stretch 




4 











CH3OH 


1128 


9.707 


35 


0.275 


- 


v-j 


CH3 rock 




4 











CH3OH 


2823 


3.543 


28 


0.035 


- 


V2 


C-H stretch 




2 











NH3 


1111" 


8.994° 


115° 


0.934 


0.6° 


V2 


umbrella 




2 











CH3OH 


1029 


9.720 


29 


0.274 


- 


Vi 


C-O stretch 




2 











CH3OH 


1132 


8.833 


35 


0.273 


- 


V7 


CH3 rock 




2 











CH3OH 


2820 


3.546 


26 


0.033 


- 


V2 


C-H stretch 




1 











NH3 


1086 


9.209 


137 


1.166 


0.8 


1^2 


umbrella 




1 











CH3OH 


1029 


9.716 


26 


0.246 


- 


Vi 


C-O stretch 




1 











CH3OH 


1135 


8.813 


44 


0.342 


- 


V, 


CH3 rock 




1 











CH3OH 


2817 


3.550 


26 


0.033 


- 


V2 


C-H stretch 




0.5 











NH3 


1080 


9.258 


118 


1.015 


0.8 


V2 


umbrella 




0.5 











CH3OH 


1030 


9.711 


22 


0.207 


- 


Vi 


C-O stretch 




0.5 











CH3OH 


1128° 


8.865° 


35° 


0.275 


- 


V, 


CH3 rock 




0.5 











CH3OH 


2813 


3.555 


27 


0.034 


- 


V2 


C-H stretch 




0.25 











NH3 


1078 


9.278 


98 


0.845 


0.9 


V2 


umbrella 




0.25 











CH3OH 


1030 


9.707 


16 


0.151 


- 


Vi 


C-O stretch 




0.25 











CH3OH 


_a 


_a 





a 


- 


V7 


CH3 rock 




0.25 











CH3OH 


2808° 


3.561° 


17° 


0.022 


_a 


V2 


C-H stretch 




1 


1 








NH3 


1116° 


8.961 


95 


0.764 


0.7 


V2 


umbrella 




1 


1 








CH3OH 


1026 


9.745 


29 


0.276 


- 


Vi 


C-O stretch 




1 


1 








CH3OH 


1125° 


8.888° 


32° 


0.253 


- 


V, 


CH3 rock 




1 


1 








CH3OH 


2824 


3.541 


26 


0.033 


- 


V2 


C-H stretch 




0.25 


10 








NH3 


1119 


8.937 


59 


0.472 


1 


V2 


umbrella 




0.25 


10 








CH3OH 


1017 


9.833 


22 


0.213 


- 


Vi 


C-O stretch 




0.25 


10 








CH3OH 


-° 


_a 


-° 


-° 


- 


V, 


CH3 rock 




0.25 


10 








CH3OH 


2829° 


3.534° 


30° 


0.037 


- 


V2 


C-H stretch 




1 


10 








NH3 


1123 


8.903 


61 


0.484 


1 


V2 


umbrella 




1 


10 








CH3OH 


1022 


9.784 


24 


0.230 


- 


Vi 


C-O stretch 




1 


10 








CH3OH 


-° 


a 


-° 


-° 


- 


V, 


CH3 rock 




1 


10 








CH3OH 


2830 


3.533 


15 


0.019 


- 


V2 


C-H stretch 




4 


10 








NH3 


1130 


8.848 


62 


0.489 


- 


V2 


umbrella 




4 


10 








CH3OH 


1023 


9.777 


30 


0.288 


- 


Vi 


C-O stretch 




4 


10 








CH3OH 


1124 


8.896 


23 


0.183 


- 


Vt 


CH3 rock 




4 


10 








CH3OH 


2830 


3.534 


14 


0.017 


- 


V2 


C-H stretch 
















CH3OH 


1028 


9.729 


28 


0.265 


1 


Vi 


C-O stretch 
















CH3OH 


1125 


8.888 


34 


0.269 


1 


V, 


CH3 rock 
















CH3OH 


2828 


3.536 


33 


0.041 


1 


V2 


C-H stretch 







1 








CH3OH 


1025 


9.755 


33 


0.314 


- 


Vi 


C-O stretch 







1 








CH3OH 


1124 


8.897 


40 


0.317 


- 


Vt 


CH3 rock 







1 








CH3OH 


2828 


3.536 


23 


0.029 


- 


V2 


C-H stretch 







9 








CH3OH 


1020 


9.801 


23 


0.221 


- 


Vi 


C-O stretch 







9 








CH3OH 


1126 


8.883 


13 


0.103 


- 


Vt 


CH3 rock 







9 








CH3OH 


2828 


3.536 


23 


0.029 


- 


V2 


C-H stretch 










9 





CH3OH 


1034 


9.675 


25 


0.229 


- 


Vi 


C-O stretch 










9 





CH3OH 


1119 


8.938 


30 


0.242 


- 


Vt 


CH3 rock 










9 





CH3OH 


2831 


3.532 


- 


- 


- 


V2 


C-H stretch 










9 





CO 


2138 


4.677 


7 


0.014 


- 


v\ 


C-O stretch 
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Table B.3 — Continued 



Ice mixture Molecule Peak position FWHM ^rci. Mode 



NH3 


CH30H 


H2O 


CO 


C02 




cm-^ 


p,m 


cm-^ 


/im 














1 
1 
1 
1 

















CH30H 
CH30H 
CH30H 

CO 


1029 
1124 
2830 
2136 


9.720 
8.898 
3.534 
4.682 


30 
32 

9 


0.286 
0.258 

0.020 


1^4 

1^1 


C-O stretch 
CH3 rock 
C-H stretch 
C-O stretch 








9 
9 
9 
9 
















CH30H 
CH30H 
CH30H 

CO 


1028 
1125 
2824 
2135 


9.730 
8.890 
3.541 
4.685 


28 
32 

9 


0.261 
0.255 

0.021 




C-O stretch 
CH3 rock 
C-H stretch 
C-O stretch 



^Band is weak and spectral overlap prohibits accurate fitting. 
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C. Comparison between astronomical and laboratory data for all sources 




2.5 3.0 3.5 4.0 4.5 5.5 6.0 6.5 7.0 

Wavelength (fLm) 



.5 9.0 8.5 10.0 




5.5 6.0 6.5 7.0 
Wavelength (^m 



0.5 10.0 



Fig. C.l. — (a) Comparison between astronomical and laboratory data for sources whose silicate absorption feature was 
fitted with a template. For a given source (displayed in either the left or right column of the figure), the middle and right 
panels show 5.2-7.5 and 8.2-10.2 fitn regions from IRS Spitzer spectra overlaid with laboratory spectra, scaled to the 9-/im 
NH3 umbrella mode, and smoothed to the Spitzer resolution. Error bars for the Spitzer spectra are indicated in the bottom- 
right corner. The dark blue line represents the pure water laboratory spectrum scaled to the water column density taken 
in paper I. Other colors are representative of laboratory spectra obtained for the following mixtures: H20:NH3=9:1 (green), 
H2 O:CH3OH:NH3=10:0 25:l (orange), H2O:CH3OH:NH3=10:l;l (cyan), and H2O:CH3OH:NH3=10:4:l (red). When available 
fsee iBoogert et al.ll200Sl) , VLT or Keck data (2.0-4.5 /im, left panel) are also plotted. In this case, we overplotted (dotted 
purple line) a pure water spectrum scaled to the 3-/im water feature of the mixed ice spectrum. Whenever present, a red dashed 
line in the right panel of a given source represents a H20:CH30H=9:1 laboratory spectrum scaled to the 9.7-/im CH3OH 
CO-stretch mode: this gives an indication of the contribution of the 9-fiBi CH3OH CHs-rock mode to the total 9-/im feature. 
The laboratory spectra are recorded at 15 K unless indicated differently. 
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A ff^ 




w 


'f.^' 




.3 3.5 




^ f 




\! 


'r 


^ 


.3 3.5 



2.0 2.5 3.0 3.5 4.0 4,5 



O.b 






1/ 




0.5 
1.0 




1 5 


.3 3.5 





VT/ 


'^ 


V 


__. 




.3 3.5 


"~^~f^ 




Y 




--»;**" 




.3 3.5 



5.5 6.0 6.5 7.0 6.5 9.0 

Wavelength ^rvi) 



a.5 3.0 3.5 4.0 4.5 




5.5 6.0 6.5 7.0 6.5 

Wavelength (/tm) 



Fig. C.l. — (b) As (a) but for sources with no associated template, i.e. with the 10-/^m sihcate feature subtracted via the 
local continuum method. Additionally, yellow represents H20:NH3=4:1 (H. Frazer, priv. comm.). 
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